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Abstract 

We study the exclusive semi-tauonic B decay, B — > Drv T , in two-Higgs-doublet 
models. Using recent experimental and theoretical results on hadronic form fac- 
tors, we estimate theoretical uncertainties in the branching ratio. As a result, we 
clarify the potential sensitivity of this mode to the charged Higgs exchange. Our 
analysis will help to probe the charged Higgs boson at present and future B factory 
experiments. 
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1 Introduction 



Many interesting models for the new physics beyond the standard model (SM) have been 
considered. One of the most attractive models is the minimal supersymmetric standard 
model (MSSM) In the MSSM, two Higgs doublets are introduced in order to cancel 
the anomaly and to give the fermions masses. The introduction of the second Higgs 
doublet inevitably means that a charged Higgs boson is in the physical spectra. So, it is 
very important to study effects of the charged Higgs boson. 

Here, we study effects of the charged Higgs boson on the exclusive semi-tauonic B 
decay, B — > Dtv t , in the MSSM. In a two-Higgs-doublet model, we have a pair of charged 
Higgs bosons, H ± , and its couplings to quarks and leptons are given by 

C H = (2v / 2G F ) 1/2 [Xu L V KM M d d R + Yu R M u V KM d L + ZU L Mil R ] H + 

+h.c. , (1) 

where M u , Ma and Mi are diagonal quark and lepton mass matrices, and Vkm is Kobayashi- 
Maskawa matrix ||. In the MSSM, we obtain 

X = Z = tan/3, Y = cotp, (2) 

where tan/5 = v%jv\ is the ratio of the vacuum expectation values of the Higgs bosons. 
Since the Yukawa couplings of the MSSM are the same as those of the so-called Model II 
of two-Higgs-doublet models ||, the above equations and the following results apply to 
the latter as well. 

From these couplings, we observe that the amplitude of charged Higgs exchange in 
B — ► Dtv t has a term proportional to m& tan 2 (3. Thus, the effect of the charged Higgs 
boson is more significant for larger tan/3. 

In Sec. 2, we give formula of the decay rate. The employed hadronic form factors are 
described in Sec. 3. In Sec. 4, we show our numerical results. Sec. 5 is devoted to conclusion. 



2 Formula of the decay rate 

Using the above Lagrangian in Eq.(|I]) and the standard charged current Lagrangian, 
we can calculate the amplitudes of charged Higgs exchange and W boson exchange in 
B -> Dtv t . 

The W boson exchange amplitude is given by [|J 

M X /(q\x) w = ^V dt Y,Kwft w m w , (3) 
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where q 2 is the invariant mass squared of the leptonic system, and x = ps • Pt/ttib- The 
r helicity and the virtual W helicity are denoted by A r = ± and \w = ±, 0, s, and 
the metric factor rj\ w is given by rj± = rjo = —r] s = 1. The hadronic amplitude which 
describes B — > D W* and the leptonic amplitude which describes W* — > tv are given by 

H{ w {q 2 ) = e;(\ w )(D(p D )\cYb\B(p B )) , (4) 
I% w {q\x) = e^(A w )(T(p r ,A T )^(p,)|T7 M (l-75)^|0> , (5) 

where e M (Av^) is the polarization vector of the virtual W boson. 
The charged Higgs exchange amplitude is given by || 



M^{q\x) H = %V cb L^ XZ*^W R + YZ*=Hl . ((>, 



Mfj K ' M 2 



v/2 

Here, the hadronic and leptonic amplitudes are defined by 



H s r ,l{<1 2 ) = (D(p D )\c(l± l5 )b\B(p B )), (7) 
L^(q\x) = (t( Pti \ t )9 t ( P „)\T(1-j 5 )v t \0) . (8) 

These amplitudes are related to the W exchange amplitudes as 

H R>L = ^—H S S , L x t = ^L^, (9) 
m b — m c m T 

where the former relation is valid in the heavy quark limit. 

Using the amplitudes of Eqs.(|3]) and (^), the differential decay rate is given by 
dV GP F \V A \ ¥V5^ 



f 1287r 3 m| 



of 



2 / q 2 tan 2 (3 m b q 2 m c ,\ 2 ( rrs^ , in , 

where Q± = (tub ± ttid) 2 — q 2 and v = a/1 — m 2 /q 2 . Note that if tan/? > 1, in which we 
are interested, this decay rate is practically a function of tan/3/M# because the second 
term in the coefficient of (H*) 2 is negligible for m^tan 2 (3 ^> m c . 

3 Hadronic form factors 

In order to obtain the decay rate numerically, it is necessary to calculate the hadronic 
amplitude in Eq.(|j). This amplitude is given in terms of hadronic form factors: 

(D(p D )\cYb\B(p B )} = ^^[h + {y){v + v'Y + h-(y)(v - v'f] , (11) 
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Figure 1: The ratios B and B as functions of R in the MSSM and the SM. The shaded 
regions show the predictions with the error in the slope parameter p\ in Eq . (|14D . The flat 
bands show the SM predictions, (a) B: the decay rate normalized to T(B — > Dpv^su- 
(b) B: the same as (a) except that the denominator is integrated over m 2 < q 2 < 
(m B - m D ) 2 . 



where v = Pb/™>b, v' = Pd/^d and y = v ■ v' — {m 2 B + m 2 D — q 2 ) /2mBm D 0. 

In the heavy quark limit and in the leading logarithmic approximation, h + (y) and 
h~(y) are given as 

h+(v) = £(v), h-(y) = 0, (12) 

where £(y) is the universal form factor 0. 

The form of £ (y) is constrained strongly by the dispersion relations as M 

£(y) ~ 1 - %p\z + (h\.p\ - 10. )z 2 - (252.p? - 84.)z 3 , (13) 

where z = (y/y + 1 — V / 2)/(a/|/ + 1 + a/2)- To determine the slope parameter pf, we use 
the experimental data of Belle 0, and we obtain 

p\ = 1.33 ± 0.22 . (14) 

This error of p\ dominantly contributes to the uncertainty in the theoretical calculation 
of the branching ratio. 



4 Numerical results 



Now, we consider the following ratio, 
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Figure 2: (a) The contour plot of upper bound of R at 90% CL as a function of Spl and 
SB exp . (b) The same as (a) except that the ratio in Eq. fllB"!) is used. 



where the denominator is the decay rate of B — > Dpv^ in the SM, since the uncertainties 
due to the form factors and other parameters tend to reduce or vanish by taking the ratio. 

Fig. 1(a) is the plot of our predictions of the ratio in Eq.(^) as a function of R, 
which is defined by R = m w ta,nj3/mH- The shaded regions show the MSSM and SM 
predictions with the error in the slope parameter pi in Eq.(|l4]). As seen in Fig. 1(a), when 
R reaches about 32, the branching ratio in the MSSM becomes the same as the one in 
the SM. It is because the interference of the W exchange and the charged Higgs exchange 
is negative. From Fig. 1(a), we expect that the experimentally possible sensitivity of R is 
~ 10, provided that the error in pi will not change. 

In Fig. 1(b), we also show the ratio, 

T(B -> Dtv t ) 



B 



T(B Dpu. 



(16) 



H)SM 



the same as Fig. 1(a), but its denominator is T(B — > Dpv^sM, which is integrated over 
the same q 2 region as the r mode, i.e., m 2 < q 2 < {tub — mo) 2 ■ From Fig.l(b), we expect 
less theoretical uncertainty and a better sensitivity of B compared with B in Fig. 1(a). 

Once the experimental values of B (B), its error SB (SB), pi and its error Spl are 
given, we can obtain a bound on R. In the following, we assume the SM prediction as the 
experimental value of B (B) , i.e., B exp = Bsm ± SB exp (B exp = B$m ± SB exp ), and we 



use the central value of Eq.flliD as the input of the slope parameter. 

Fig. 2(a) is the contour plot of upper bound of R at 90% CL as a function of Spl 
and SB exp . From this figure, if SB exp = 0, and Spl ~ 17%, which corresponds to the 
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present experimental error in Eg. flT4]) , we expect that an upper bound of R ~ 10, which 
is consistent with the result of Fig. 1(a). If we will observe B — > Dtv t with 5B exp ~ 20%, 
we expect that an upper bound of R ~ 15 weakly depending on 5p\. 

In Fig.2(b), we also show a similar contour plot where we use the ratio defined in 
Eq.(|l6|), i.e., normalized to T(B — > D^u^sm- We observe that the upper bound of R is 
almost independent of 5p 2 in this case. Thus, it is important to make the experimental 
error in B (B), 5B exp (5B exp ), small rather than 5p 2 . 



5 Conclusion 



As seen in our numerical results, the branching ratio of B — > Dtv t is a sensitive probe of 
the MSSM-like Higgs sector. We expect an upper bound of R < 15 when 6B exp ~ 20% 
is achieved. So, if B — > Dtv t is observed at a B factory experiment, a significant regions 
of the parameter space of the MSSM Higgs sector will be covered. Comparing with the 



Higgs search scenario of LHC |10|, we conclude that present and future B factories are 
potentially competitive with LHC. 

As future improvements of the present work, the q 2 distribution [] 1] and the r polar- 
ization || of B — > Dtv t are promising. In these quantities, we will expect that the theo- 
retical uncertainties from the error in the slope parameter become very small. However, 
we should take 1/m and QCD corrections into account. These corrections are neglected 
in the present work because they lead to smaller uncertainties than those from 5p\. For 
the q 2 distribution and the r polarization, they are expected to be dominant uncertainties 
in the theoretical calculations. These issues will be addressed elsewhere. 
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